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ABSTRACT 


The Philippine sardines (Genus Sardinella, Family Clupeidae) are interest- 
ing organisms due to their great morpho-anatomical similarities. They usually 
thrive in marine environment with the exception of the freshwater sardine 5. fawilis: 
endemic to Lake Taal. To elucidate the relationship of the Philippine Sardinella 
‘pp, 35 biometric variables were determined and then subjected to descriptive 
(mean, standard deviation, and statistical range) and inferential (Factor Analysis 
(FA), Discriminant Function Analysis (DFA), and Hierarchical Cluster Analysis 
(НСА) tests. Thirty variables correlated to the four-factor data reduction in FA 
that explains 96.9% of the total variance among the six sardines. In DFA, only 29 
variables were found useful in differentiating the sardines. The same 29 features 
мете able to properly predict the identity of 99.58% of the sardines. The close 
relationship of S: albe/la to the Lake Taal sardine was demonstrated by all three 
dendrograms drawn from all 35 variables, only from the 30 FA informative vari- 
ables, and from the 29 discriminating variables, respectively. The dendrograms 
also show that the two populations of S. awilis form two distinct clusters consis- 
lent with evidence from molecular data. Among the marine sardines, S. lemuru 
сап be regarded as the farthest relative of the freshwater sardine S. sawilis. 
Keywords: Sardinella, Clupeidae, morphometrics, meristics, descriptive and 

multivariate tests 





INTRODUCTION 


Sardines are moderately small, silvery fishes that live in the shallow depths 
of coastal water or occasionally in estuaries. They are known for their oily flesh 
‘and their two-lobed non-functional lung (Pileggi and Thompson, 1979). In addi- 
tion, sardines have thin, deciduous scales that cover their elongate, spindle shaped 
bodies. These fishes comprise the genus Sardinella that belongs to the family 
Clupeidae [Order Clupeiformes; Class Actinopterygii]. They are considered one 
ıı the most important food fishes in the world, providing food not only to man, 
but also to large marine predators and sea birds (Herald, 1961). In the Philip- 
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pines, there are seven species identified under genus Sardinella, six of which are 
marine species (S. albella, Saal: 5. brachysoma, Sabr; S. fimbriata, Safi; 5. lemuru, 
Sale; S. longiceps, Salo, and S. melanura) and the other one is a freshwater 
species (S. tawilis, Saal) that can only be found in Taal Lake (Conlu, 1986). 

S. tawilis was formerly named Harengula tawilis but was now given its new 
name based on its morpho-anatomical similarity with the members of the Sardine lla 
genus. Known to local folks as tawilis, it was given an English name of freshwa- 
ter herring. Tawilis can be caught throughout the year particularly during the 
months of September and November using gill nets set with a float-line found at 
а short distance below the surface of the lake (Herre, 1927). Other members of 
the genus like $ a/bella (white sardinella), S. brachysoma (deepbody sardinella), 
S. fimbriata (feingescale sardine), 5. lemuru (Bali sardinella), 5. longiceps (In- 
dian oil sardine), and 5. melanura (Blacktip sardineila) are abundant in the Balayan 
Bay of the South China Sea and other bays in the country (Rau and Rau, 1980; 
Conlu, 1986). Unlike tawilis, all the other sardines are pelagic species inhabiting 
coastal waters, These are called by different names depending upon the region 
where they are being caught. To the locals, most of these sardines are called | 
Tamban. Other names are also used like Tuloy and Bagabas (Tagalog region), 
Alupaypay (Bicol), and Chevy (Sorsogon). 

Delineation of species within the genus has so far been based on morpho- 
anatomical features alone, which in many cases, is quite confusing because of 
their great similarity with each other. Taxonomists have tried to resolve this prob- | 
lem by using rigid criteria for discriminating species that can be derived through | 
descriptive, multivariate, and cluster analyses (Gnanadesikan, 1977), These meth- | 
ods have been already been applied to facilitate fish research, Beacham (1985)! 
succeeded in differentiating the brood lines of pink salmon using meristic and 
‘morphometric characters, The same method was used by MacCrimmon and, 
Clayton (1985) to detect the regional home-river origin of salmon, In a trawl 
survey of 17 groundfish species in the Ѕошћет Adriatic Sea, biometrics 
become useful in elucidating a new species of Caelorinchus (Ungaro et al., 2001); 
‘The phylogenetic relationship of the diverse and complex genus Barbus was re 
fined using biometric features of the bones, cleithra, opercular, and pharyngeal 
(Escala and Miranda, 2003). Through biometrics, populations of Cobitis 
‘simplicispina and C. turcica from the Sakarya-Kizilirmak river basins, Bagilli 
and Saz Lake in Turkey were were established as unique populations (Erk’ahan 
е! al., 2003). Sixteen morphometric and eight meristic characters, on the other 
hand, contributed significantly in the differentiation of the two populations of 
Sabanejewia balcanica inhabiting two rivers in Central Croatia (Bucar ef al, 
2003), In sardines, Intrina (1999) has employed the use of biometrics to describe 
three (Saal, Safi, and Sata) of the seven Philippine species. 

‘This study was performed with the general aim of resolving the relationship 
among the different Philippine sardines, with the exception of S. melanura that is 
difficult to find and identify, using a set of biometric characters/features, 
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MATERIALS AND METHODS 


Collection and Identification of Fish Samples 

Five Philippine sardines [5. albella (Saal), 5, brachysoma (Sabr), $. lemuru 
(Sale), 5. longiceps (Salo), and S. tawilis (Saraj) were collected and transported 
1o the Science and Technology Research Center (STRC), De La Salle University 
101.50) – Manila. The first four species were bought from the Navotas Fish Port, 
Navotas, Metro Manila. Interview with fish dealers revealed that the fish was 
гот Cebu City and the neighboring Batangas Вау. Two populations of S. tawilis 
(Sata) designated as “SataA” and “Sara” were collected from the northern- 
(Talisay) and the southern- (Agoncillo/San Nicholas) basins, respectively, of 
Lake Taal. 

‘The identity of the sardines was ascertained using the Key to the Species of 
Clupeinae in the Westem Pacific (Munroe ef al, 1999, Fishbase 2005). 


Thiny-five (35) individuals each of pre-identified Saal, Sabr, Sale, and the 
бо populations of Sara were used for biometrics. Only 26 individuals of Salo, 
оп the other hand, were tested due to low catch of this species during the time of 
sampling, 

Biometric Methods 


A total of thirty five variables consisting of three meristic (numbers of dor- 
sals, pectoral-, and ventral fin-rays) and 32 morphometric variables (Table 1) 
were counted and measured in the sardines. 

Measurements were made using Vernier caliper (Mitutuyo) calibrated to the 
‘nearest 0.01 mm. In cases wherein measurements exceeded the capacity of the 
caliper, a millimeter ruler was used. Bilaterally symmetrical features were mea- 
sured or counted on the morphometric landmarks on the left side of the fish as 
shown in Figure 2 and defined in Table 2. Nomenclature for the inter-landmark 
stances (Table 1) followed that of Armbruster (2003). 

Data Analyses 
Meristic and Morphometric 


Biometrics data from the five Philippine sardines were combined with the 
ХИЛ data of Intrina (1999), The raw data was obtained from 236 individuals be- 








let measures of central tendency (mean, standard deviation, statistical range), 





variate tests (factor- and discriminant function analyses), and hierarchical 
‘luster analysis, 

‘The measures of central tendency (mean, standard deviation, and statistical 
те) were computed using the statistical functions in Microsoft Excel, The lower 
Мон (LL) of the statistical range is computed as the difference between the mean 
ıl twice its standard deviation (LL = mean - 2SD), The upper limit (UL), on the 
‘ulir hand, is equal to the sum of the mean and twice its Standard deviation (UL 
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= mean + 250). The ranges of the meristic measurements were rounded off to the! 
nearest whole number. 

For multivariate analyses, STATISTICA (version 7.0) was used. In Factor 
Analysis (FA), eigenvalues and their corresponding normalized eigenvectors were 
derived from the correlation matrix of the data. These eigenvalue-cigenvector 
pairs served as the principal components (PCs) in the analysis. The principal 
components whose eigenvalues were greater than or equal to one (E е” 1) were 
retained as significant factors, according to Kaiser's Rule The factor pattern was 
rotated using the varimax raw method to magnify the proximity of a single vari- 
able with a particular factor while minimizing its proximity tothe other factors, A 

jon value/factor loading of е” 0.700 is needed for a variable to be in- 





Discriminant Function Analysis (DFA), on the other hand, was used to de- 
sign a model that discriminates the different species following the forward stepwise 
method. The variables were added one by one in the model. The ability of the 
variables to differentiate the specified groups is represented by F values. An F 
value is che iu © the between-groups variance in the data over the average 
satnin-group vce. If the between-group variance is significantly larger than 
the variance within the group then there must be significant differences between 
the means, In this test, measurements with an F value greater than or equal to 1 (Е 
e” 1) contribute significantly in discriminating one species from the others. How- 
ever, those variables with an F value of less than 1 (Е < 1) do not. The constructed 
model was then allowed to run through the data to predict each sample to its 
corresponding species. 

Hierarchical Cluster Analysis (HCA) was performed in PCord (version 3.20) 
using Centroid method as the linkage rule. Centroid method is an agglomerative 
algorithm in which the similarity is measured as the distance between cluster 
centroids, The centroid of a cluster is the average value of the objects obtained in 
the cluster on each variable. When two clusters are combined, a new centroid is 
computed. Due to the enormous data available, only the descriptive mean values 
of the sardines were utilized in establishing the relationship of the six sardines. 
‘Three dendrograms, therefore, were constructed that made use of data from: afl 
35 variables, only the correlated variables in FA, and only the discrimi 
variables in DFA. і 





RESULTS AND DISCUSSION 
Descriptive Statistics on Sardines’ Biometric Data 

‘Measurements of three meristic and 32 morphometric variables from а total | 
of 236 individuals belonging to 6 different species of Philippine sardines com- | 
prised the entire raw data. With the exception of Safi data (kindly provided by 
Ms. Marnie Grace Intrina-Sonico), the rest of ће measurements for the 2 popula- 
tions of Sata, Saal, Sabr, Sale, and Salo are generated in the present study, 


36 


\ulippine Journal of Systematic Biology Vol i, No. 1 (June 2007) 





The mean, standard deviation, and statistical range were computed in all the 
meristic (Tables 3 to 5) and morphometric data (not shown). Table 3 presents а 
comparison of the sardines’ dorsal fin ray counts obtained from different studies 
In the present study, the dorsal fin rays of sardines ranges from 12 to 19. This 

alue somehow fits im the 13 to 21 statistical range previously reported (White- 
head, 1985), 

‘The pectoral fin rays of sardines range from $ to 17 (Table 4). S. albella and 
У. lemuru have 14 to 16 fin rays, whereas, the counts for 5. sawilis, $. brachysoma, 
ind S. fimbriata range from 8 to 15 Due to the relative wide range pectoral fin 
тау counts in sardines, the number of pectoral fin rays may not be a good criterion 
Тог distinguishing the different sardines. Whitehead (1985) also did not include 
pectoral fin rays as а variable in his study. 


Another meristic character/variable considered for sardine identification was 
the ventral fin ray count. The mean with standard deviation, as well as, the statis- 
\ical ranges of the ventral fin rays of sardines are presented in Table 5 Whitehead 
(1985) reported that sardines have either 8 (S. albella, 5. brachysoma, S. fimbriata, 
and S. tawilis) ог 9 (S. lemuru) ventral fin rays. Results on 5. brachysoma (8 
ventral fin rays) and S. lemuru (9 ventral fin rays) further confirmed Whitehead’s 
findings (1985) suggesting that the number of ventral fin ray can be used to dis- 
Linguish S. lemuru from the rest of the sardines studied, The data of Intrina (1999) 
on the ventral fin ray counts of S. sawilis, 5 albella, and 5. fimbriata, on the other 
hand, yielded a range of 6 to 9 ventral fin rays 


Overall, the descriptive statistics for the 32 morphometric variables were not 
that informative. From these values alone, it is quite difficult to ascertain which 
morphometric variables can actually be useful in distinguishing the different sa 
dines. Hence, multivariate tests were conducted 

Factor- and Discriminant Function Analyses on Sardines’ Biometric Data 


In FA, four factors were found to satisfy Kaiser's rule (eigenvalue of е” 1) 
which means that the 35 variables can now be explained in terms of four smaller 
groups of variables, otherwise called as factors These factors were retained as 
the principal components (PCs) in the Principal Component Analysis (PCA) which 
account for a total of 96.9% of the total variance observed among the six sar- 
dines. 


Rotation of variables using varimax method yielded the variables correlated 
to each of the four factors previously identified. A variable must have correlation 
value/factor loading of >0.700 to be correlated to a given factor. A complete list- 
ing of variables correlated to Factors 1, 2, 3, and 4 is given in Table 6. Seventeen, 
eight, four, and one, variables, respectively, correlated to Factors 1, 2, 3, and 4 
‘The remaining five variables (P, hDe, ocDH, mL, and Lep) were found to have no 
substantial contribution to the total variance present among the samples. There- 
fore, these variables were not correlated to any of the four factors. Since all vari- 
ables under each factor are highly correlated, it is therefore possible to consider 
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‘only one variable under each factor and still give a similar variance when all 
variables were taken together 

In DFA, 29 variables were found useful in differenting the five sardines (Table 
7) All these variables have Е values greater than or equal to 1 (Е е” 1), The rest of 
the variables (рар, v-aD, V, dfL, Lvf, Lpf) have F values less than 1 (Е < 1), thus, 
were excluded in the analysis. The 29 discriminating variables were then used in 
the forward 29-stepwise Discriminant model for sardines. Each step in the modet 
is characterized by the addition of the variable with the highest contribution in 
discriminating between the specified groupings (in this case, by species). 

As shown in Table 8, the model predicted the identity to the species level of 
99.58% (235 out of 236 individuals) sardines, The model works perfectly (100%) 
in classifying Sata, Sabr, Safi, Sale, and Salo. The only misclassification was 
observed in Saall05 which the model identified as SataT. 





Relationships of Sardines Revealed by Biometrics Data 


The three dendrograms constructed based on (1) all 35 biometric variables, 
(2) 30 highly correlated variables, and (3) 29 discriminating variables only gave 
similar topology and is therefore presented as one consensus tree (Figure 3). The 
sardines are distributed into five major branches [(SaraA, SataT, and Saat), (Sabr), 
(Safi), (Sale), and (Salo)]. What is interesting in this clustering is that S. abella 
(Saal) grouped with the two populations of Sata. This grouping supports the pre- 
vious findings based on molecular data that S. a/bella is the most likely ancestor 
to the Taal Lake sardines (Samonte ег al, 2000) with S. femuru as the farthest 
relative. Also, the separation of SalaA and SataT populations into two distinc 
‘groups indicates genetic differentiation that at this point can already be observ 
morphologically. The differences may have developed in each of the popula 

due to reduced fish migration across a shallow barrier in the middle of the 
Conclusion 








‘The relationship of Philippine sardines was determined using biometrics. A 
total of 35 biometric variables (3 meristic and 32 morphometric) were determined 
оп 236 individuals of sardines belonging to six different species. The biometric 
data were then subjected to descriptive (mean, standard deviation, and statistical 
tange) and inferential (FA, DFA, and HCA) tests. 

From the descriptive values obtained, the ventral fin ray count appears to 
а promising feature that distinguishes S. /emuru from the rest of the sardi 
congruent to earlier findings of Whitehead (1985), 

1n FA, the 35 variables are reduced to 30 variables only that are highly corre- 
lated to four factors (eigenvalue > 1). These factors or small groups of variables 
account for a total of 96.9% of the total variance observed among the six 
dines. Twenty nine (29) variables, on the other hand, were found useful in dis- 
ctiminating the Philippine sardines. These variables were used in the 29-stepwise 
discriminant model that correctly predicted identity of 99.58% (235 out of 236) 
individual sardines. 
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CONCLUSIONS 


Finally, hierarchical cluster analysis (based on all 35 biometric variables, 
only the 30 highly correlated variables, and the 29 discriminating variables) ге- 
vealed that phenetically, S. albella is the closest marine relative of the freshwater 
У tawilis with S. lemuru as the most distant marine associate, The separation of 
the two populations of $. tawilis (SataA and SataT) is also consistent with earlier 
findings based on molecular data (Samonte ef a, 2000). Thus, supporting the 
hypothesis of an on-going genetic differentiation between the two populations. 
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Figure 1. Sardinella spp used in this study. 
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Figure 2, Landmarks and distances measured. Broken lines indicate per- 
pendicular lines transecting the total body length (1-27) from their corre- 
sponding point of reference. See Table 1 for the 











morphometrics and Table 2 for definitions of morphometric landmarks. 
Constructed based on Armbruster, 2003. (Drawn by К. Alvia). 
Distance (Objective Function) 
БЕД REXI 2020,00 4226.429 5624.429 
Пн چن‎ — = 
Information remaining (0) 
100.000 75.000 50.000 25.000 000 





Sale 





Figure 3. Representative dendrogram showing the relationship of the 
different sardines using Centroid linkage rule method. 
a 


Philippine Journal of Systematic Biology Vol. i, Хо 1 (June 2007) 





‘Table 1. Morphometric variables used in this study. Landmarks repre- 











sent the two points from where measurements were made in Fig. 2. 
Enclosed in parentheses are abbreviations of the variables. 
Landmark Variable 
12 pre orbital length (preOL) 
23 horizontal diameter of the eye (De) 
34 post orbital length (posOL) 
14 head length (hL) 
56 orbital depth of the head (orDH) 
78 occipital depth of the head (ocDH) 
9-11 lower jaw length (jL) 
10-11 maxillary length (mL) 
12-13 maxilliary width (mW) 
14-15 dorsal fin base length (dfL.) 
14-16 length of the first ray of the dorsal fin (Lladf) 
17-18 length of the fourth ray from the posterior end of the dorsal 
fin (L4pdt) ү 
19-20 length of the third ray from the posterior end of the dorsal fin 
(02540 
21-22 length of the upper lobe of the caudal fin (Lulcf) 
23-24 length of the middle lobe of the caudal fin (Lmlct) 
25-26 length of the lower lobe of the caudal fin (Lilet) 
1-27 ‘maximum total length (tL) 
1-24 fork length (fL) 
1-28 ‘maximum standard length (sL) 
14-29 maximum body depth (maxbd) 
30-31 ‘minimum body depth (minbd) 
1-33 preanal distance (paD) 
1-35 preventral distance (pvD) 
136 predorsal distance (рар) 
35-37 pectoral-ventral distance (p-vD) 
35-40 ventral-anal distance (v-aD) 
28-40 anal-caudal distance (aD) 
28-44 length of caudal peduncle (Lep) 
4142 height of anal fin (Haf) 
3943 length of anal fin base (Lafb) 
3445 length of ventral fin 1м) 
3746 length of pectoral бп (Lpf) 
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ble 2, Definitions of morphometric landmarks in Figure 2. 





Landmark Definiti 








T  tipofmouth 

2 anterior border of the left eye 

3 posterior border of the left eye 

4 edge of operculum 

5 dorsal margin of the orbital depth of the head 

6 ventral margin of the orbital depth of the head 

7 dorsal margin of the occipital depth of the head 

3 ventral margin of the occipital depth of the head 

9 tipof lower jaw 

10 tip of upper jaw 

11 junction of the upper and lower jaws 

12 right border of upper jaw 

13 left border of upper jaw 

14 anterior end of dorsal fin base 

15 posterior end of dorsal fin base 

16 tip of the first dorsal fin ray 

17 insertion of the fourth ray from the posterior end of the dorsal 
fin base 

18 tip of the fourth ray from the posterior end of the dorsal fin base 

19 insertion of the third ray from the posterior end of the dorsal fin 


base 
20 tip of the third ray from the posterior,end of the dorsal fin base 
21 insertion of the upper lobe of the caudal fin 


22 tipof the upper lobe of the caudal fin 

23 insertion of the middle lobe of the caudal fin 

24 tip of the middle lobe of the caudal fin 

25 insertion of the lower lobe of the caudal fin 

26 tip of the lower lobe of the caudal fin 

27 point in line / perpendicularly transecting the tip of the lower 
lobe of the caudal fin 

28 point in line /t perpendicularly transecting the insertion of the 
lower lobe of the caudal fin 

29 ventral margin of the maximum body depth 

30 dorsal margin of the minimum body depth 

31 ventral margin of the minimum body depth - 

32 anal opening 
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‘Table 2. (continued) 1 





34 
35 


36 


37 
38 


39 
40 


s 
42 
43 
44 


45 
46 


point in line /t perpendicularly transecting the anal opening 
anterior end of the ventral fin base 

point in line /¢ perpendicularly transecting the anterior end of the 
ventral fin base 

point in line И perpendicularly transecting the anterior end of the 
dorsal fin base 


anterior end of the pectoral fin base 


point in line /t perpendicularly transecting to the anterior end of 
the pectoral fin base 


anterior end of the anal fin base 


point in the line / perpendicularly transecting the anterior end of 
the anal fin base 


proximal end of the anal fin height 
distal end of the anal fin height 
posterior end of the anal fin base 


point in the line / perpendicularly transecting to the posterior end 
of the anal fin base } 


tip of the first ray of the ventral fin 
tip of the first ray of the pectoral fin 











Note: 1, total body length 


Table 3. Descriptive values ов the dorsal fin ray count of fivellpecies of sar. 
dines, Ranges in the present study are computed as mean + 2SD roundea 
off to the nearest whole number. 

















CESS ee 4 
Sarge n Мн? O Preet зв Whitehead 
study _ 01999) (1955) 
Samak 55 1732072 } 
жшт иа jer (ра 
Seat 35 158142 15-18 13-21 
Sabr 35 1694048 4 13-11 
Saft 35 172401 16-18 43-21 
Sale 35 1592049 д 
1542058 














не а mumber of individuals per species. SD. standard deviation 
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Fable 4. Descriptive values on the pectoral fin ray counts of five species of 
sardines. Ranges in the present study are computed as mean + 2SD 
rounded off to the nearest whole number. 





Range 








Sardine ХО Mean=SD Present Inivina 
study _ (1999) 

Saran 14321 1 13-15 

Sarat 12 s - 

Sani Ш7=173 8-15 14-16 

Sadr 141=037 13-15 

5ай 1452051 14-16 13 

Sue 1S3=048 14-16 

Saio 26___143073__1316 





Note п. шипбет of mdividnals per species SD standard deviation 


‘Table 5. Descriptive values on the ventral fin ray counts of five species of 
sardines, Ranges in the present study are computed as mean + 2SD 
rounded off to the nearest whole number. 














Range 
Sardine n Meam=SD рн  Intrina Whitehead 
study (1999) (1985) 

Saran 1 

sî 3 r + 
Sal 35 7.9 3 
Sar 35 8 
Safi 55 79 8 
Saie 35 5 э 
Salo 26 769=047 1-9 3 





Nine п. umber of individuals per species. SD, Standard deviation 


45 
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‘Table 6. Summary results of Factor Analysis showing reduction to 
four factors and their corresponding correlated variables 





Factor Correlated Variables 





1 V, hL, preOL, posOL, orDH, mW, L4pdf, L3pdf, tL, fL, sL, рар, 
`pvD, рар, p-vD, v-aD, Lpf 


2 diL, Lulef, Lilef, maxbd, minbd, a-cD, Lafb, Lvf 
3 D, jL, тіс, Haf 
لیے‎ 





‘Table 7, List of discriminating and non-discriminating variables 











Variables 
Discriminating Гара, iL, a-cD, minbd. Lmicf, Lladt, Haf, 
Variables posOL, maxbd, D, mW, mL, Lafb, Lep, p-vD, 1 


AL, orDH, fL, Lulcf, hL, рар, осн, pvD,hDe, 
3pdf, рар, preOL, P.Lilef, SL 
Non-discriminating 
Variables рар, v-aD, V, dfL, Lvf, Lpf 





‘Table 8. Discriminant Function Analysis predictions (%) for the correct iden- 
tification of Philippine sardines. 











Observed п Sumber of dividuals predivied eon 

sardine Seiad Sara Saai Sobr Safi Sale Soio уы 
шм з з о 0 о о о 0 100.00 
заат 33 о 3 о 9 0 0 0 100.90 
Sai 5 о а ын ооо о 9714 
Sor з о о е 3 о 0 о 100.90 
Sf 035. 19. '@- OS бе 20 100.00 
Sie 5 o о о о о 35 о 100.00 
Se 26 100 00 
Toul _ 236 эз 





‘Note * incorect DFA-model identifications 


